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BAYESIAN APPROACH TO INVERSE PROBLEMS FOR 
FUNCTIONS WITH VARIABLE INDEX BESOV PRIOR 


JUNXIONG JIA, JIGEN PENG, AND JINGHUAI GAO 


Abstract. We adopt Bayesian approach to consider the inverse problem of 
estimate a function from noisy observations. One important component of 
this approach is the prior measure. Total variation prior has been proved 
with no discretization invariant property, so Besov prior has been proposed 
recently. Different prior measures usually connect to different regularization 
terms. Variable index TV, variable index Besov regularization terms have 
been proposed in image analysis, however, there are no such prior measure in 
Bayesian theory. So in this paper, we propose a variable index Besov prior 
measure which is a Non-Guassian measure. Based on the variable index Besov 
prior measure, we build the Bayesian inverse theory. Then applying our theory 
to integer and fractional order backward diffusion problems. Although there 
are many researches about fractional order backward diffusion problems, we 
firstly apply Bayesian inverse theory to this problem which provide an oppor¬ 
tunity to quantify the uncertainties for this problem. 


1. Introduction 


Partial differential equations are powerful tools for describing physical systems. 
Using partial differential equations we can predict the outcome of some measure¬ 
ments. The inverse problem consists of using the actual result of some measure¬ 
ments to infer the values of the parameters that characterize the system. Partial 
differential equations usually have a unique solution, while the inverse problem 
does not [T]. Because of this, we usually need to use some a priori information to 
compensate the data uncertainties. 

Regularization techniques are useful tools to produce a reasonable estimate of 
quantities of interest based on the data available. Studies about regularization tech¬ 
niques have a long history dating back to A. N. Tikhonov in 1963 [2]. When study¬ 
ing backward diffusion problems, the researchers found that the classical Tikhonov 
regularization with standard penalty terms such as || Vu ||^2 is well known to be al¬ 
ways over smoothing. So total variation (TV) regularization has been proposed in 
[3]. TV penalty have gained increasing popularity for it can preserve the important 
details such as edges of the image. In 1997, P. Blomgren, T. F. Chan, P. Mulet 
and C. K. Wong [4] noticed that TV restoration typically exhibit “blockiness”, or 
a “staircasing” effect where the restored image comprises of piecewise flat regions. 
So they proposed a regularization term as follows 



where p monotonically decreasing from 2, when |Vm| = 0, to 1, as |Vu| oo. Later 
in 2006, Y. Chen, S. Levine and M. Rao [5] proposed another kind of variable index 
TV norm which has many mathematical structures. In 2014, J. Tiirola [6] used 


1 


2 


J.X.JIA, J. PENG, AND J. GAO 


variable index TV norm and variable index Besov regularization terms in image 
decomposition problems. 

However, regularization techniques can not give uncertainty analysis. Statistical 
inversion theory reformulate inverse problems as problems of statistical inference 
by means of Bayesian statistics. Dating back to 1970, Franklin [7] formulating 
PDE inverse problems in terms of Bayes’ formula on some Hilbert space X. In 
this paper Franklin deriving a regularization using Baysian approach, it also state 
the relation between regularization techniques and Bayesian approach. Recently, 
S. Lasanen [sjiaiiniiii] developed fully nonlinear theory. S. L. Cotter, M. Dashti. 
J. C. Robinson, A. M. Stuart, K. J. H. Law and J. Voss [27l [28l [32] establish a 
mathematical framework for a range of inverse problems for functions given noisy 
observations. They establish the relation between regularization techniques and 
Bayesian framework, estimate the error of the finite dimensional approximate solu¬ 
tions. Based on this framework, S. L. Cotter, G. O. Roberts, A. M. Stuart and D. 
White [12] develop faster MCMC algorithms. 

Now, we need to state the relation between TV regularization and Bayesian 
approach. In [13], M. Lassas and S. Siltanen find TV regularization is not dis¬ 
cretization invariant, that is to say, Bayesian conditional mean estimates for total 
variation prior distribution are not edge-preserving with very fine discretizations 
of the model space. In order to overcome this deficiency, M. Lassas, E. Saksman 
and S. Siltanen |22j proposed Besov prior B\ ^ which is discretization invariant. 
M. Dashti, S. Harris and A. M. Stuart [^ studied the Besov prior under the 
mathematical framework established in m- Under this framework Besov prior is 
discretization invariant naturally for the framework is originally build on infinite 
dimensional space. 

Considering the Besov and TV regularization techniques and Bayes’ inversion 
theory, we find that there may no variable index Besov prior theory. As mentioned 
before variable index TV and variable index Besov regularization terms have been 
used in image analysis and achieved good performance. In this paper, we try to 
build variable index Besov prior and generalize the Bayes’ inversion theory with 
this new prior probability measure. In section jb] we use our theory to integer order 
backward diffusion problems and fractional order backward diffusion problems. 

The main contributions of this paper are 

(1) Construct a variable index Besov prior using wavelet characterization of 
variable index Besov space and prove a Fernique-like result [14] for variable 
index Besov prior; 

(2) Based on variable index Besov prior, we generalize results in [28] to build the 
Bayes’ inversion theory. Under same conditions for the forward operator, 
we also prove the convergence of variational problems with variable index 
Besov regularization term; 

(3) Although there are many studies about inverse problems for fractional dif¬ 
fusion equations [dO] [41] , until now, there are few studies about fractional 
order backward diffusion problems under Bayes’ inversion framework. Us¬ 
ing our theory, we proved the posterior measure exist and the continuity of 
the posterior measure with respect to the data for integer order backward 
diffusion problems and fractional order backward diffusion problems. 

The content of this paper are organized as follows. In section [21 we state some 
basic knowledge about variable index space and prove wavelet characterization of 
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variable index Besov space with periodic domain. In section [3l we construct the 
variable index Besov prior and proved a Fernique-like theorem. In section [4l we 
generalize Bayesian inversion theory to our variable index Besov prior setting. In 
section [SJ under same conditions in section |4] for forward problem, we proved the 
variational problem with variable index Besov regularization term converge. In 
section we used our theory to integer order backward diffusion problem and 
fractional order backward diffusion problem. In the last two sections, we give some 
technical lemmas and for the reader’s convenience we list some useful theorems and 
lemmas used in our paper. 

2. Variable order space and wavelet characterization 

In this section, we give a short introduction to space of variable smoothness 
and integrability on periodic domain, then we prove a wavelet characterization of 
variable index Besov space on periodic domain. 

2.1. Modular spaces. 

Definition 2.1. [20] Let V be a vector space over M or C. A function p : A 
[0, oo] is called a semiinodular on X if the following properties hold: 

(1) p(0) = 0. 

(2) p{\f) = p{f) for all / e A and |A| = 1. 

(3) p(A/) = 0 for all A > 0 implies / = 0. 

(4) A i-A p{Xf) is left continuous on [0,oo) for every / G A. 

A semimodular p is called a modular if 

(5) Pif) = 0 implies / = 0. 

A semimodular p is called continuous if 

(6) for ever / G A the mapping A M- p{\f) is continuous on [0, oo). 

A semimodular p can be additionally qualified by the term convex. This means, 
as usual, that 


p[ef + (i-e)g)<ep{f) + {i-e)p[gi 


for all f,g G X. 

Once we have a semimodular in place, we obtain a normed space in a standard 
way: 

Definition 2.2. |20| If p is a (semi)modular on A, then 


Ap := {a; G A : 3 A > 0, p{Xx) < oo} 


is called a (semi)modular space. 

Theorem 2.3. |20j Let p be a convex semimodular on X. Then Ap is a normed 
space with the Luxemburg norm given by 
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2.2. Spaces of variable integrability. The variable exponents that we consider 
are always measurable function on n-dimensional torus T" with range [l,oo). We 
denote the set of such functions by V. We denote = esssup^g.].„p(a;) and p~ = 
essinf 2 ,gT'‘P(a^)- The function pp is defined as follows: 

( tP, ifp€(0,oo), 

‘Ppi't) = A 0, if p = oo and t < 1, 

[ oo, if p = oo and t > 1. 

The convention 1°° = 0 is adopted in order that pp is left continuous. In what 
follows we write tP instead of <Pp{t). The variable exponent modular is defined by 

PLr()(/):=/ \f{x)\^^^^dx. 

The variable exponent Lebesgue space LP^''> and its norm ||/||p(.) are defined by the 
modular as explained in the previous subsection. 

We say that p : T” ^ R is locally log-Holder continuous, abbreviated g € 
if there exists c > 0 such that 


\9i.x)-g[y)\ < 


log(e + l/|x-p|) 


for all x,y G T”. We say that g is globally log-Holder continuous, abbreviated 
g £ C’‘°P, if it is locally log-Holder continuous and there exists poo G R such that 


Iff (a:) -ffool < 


log(e-b |a;|) 


for all a; G T”. The notation is used for those variable exponents p £ V with 
1 g (jiog that is to say 1 < p“ < p{x) < p+ < oo and ^ is globally log-Holder 
continuous. 


2.3. Variable index Besov space. Before we introduce variable index Besov 
space, we need the following definition of mixed Lebesgue-sequence space. 

Definition 2.4. [17] Let p,g £V. The mixed Lebesgue-sequence space i'^^'\LP^''>) 
is defined on sequences of L*’*^')-functions by modular 

:=^inf jA„ > 0 | ) (/„/A?^) < l|. 

As usual, denote the Fourier transform of a distribution or a function / as J^{f) 
or /. Denote the inverse Fourier transform of a distribution or a function / as 
/^. As in the constant index case, we need the following definition of admissible 
functions. 

Definition 2.5. m We say a pair {ip, 4)) is admissible if $ G 5 satisfy 

• suppp C G T” : 1/2 < 1^1 < 2} and |:p(^)| > c > 0 when 3/5 < |^| < 

• supp4> C G T" : 1^1 < 2} and |$(OI > c > 0 when |^| < 5/3. 

We set (pv{x) := 2™ip{2'"x) for f G N and (po{x) := <l’(a;). Denote S to be the 
Schwartz function space. S' to be the tempered distribution that is the dual space 
of S. Then the variable index Besov space in our setting can be defined as follows. 
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Definition 2.6. [17] Let (pv be as in Definition 12.51 For a : T" —>■ R and p,q gV, 
the variable index Besov space consists of all distributions f G S' such 

that 


3“(') 




li(-) {LPi )) 


< oo. 


In the case of p = g we use the notation := ^p{-) q{-)- Besov space 

we can also associate the following modular: 

) := Pp!( )(Lp( ))((2”“^'V« * /)«)> 

which can be used to define the norm. For the reader’s convenience, we also list 
the definition of variable index Triebel-Lizorkin space ^pp)qp)- 


Definition 2.7. [53| Let tpv, i; G N U 0, be as in Definition 12.51 The Triebel- 
Lizorkin space Fp^.)q(^.') is defined to be the space of all distributions f £ S' with 
ll/ll c(.) < oo, where 

P(-).9(-) 




.(■) 


T’p * / 


LP() 


In the case of p = 
this paper, we denote 


q we use the notation := In the following of 

A K, B equal to cA < i? < CA with c, C be two constants. 


2.4. Wavelet characterization. Now, we state some notations for wavelet theory 
then prove a wavelet characterization of variable index Besov and Triebel-Lizorkin 
space on periodic domain. 

Let '0'^, 0^ be the Meyer or Daubechies wavelets described in Proposition 19.11 
in the Appendix. Now define 

G° = {F, M}" and = {F, M}’"* if j > 1, 


where the * indicates that at least one Gi of G = (Gi, • ■ • , Gn) S {F, M}"* must 
be an M. It is clear from the definition that the cardinal number of {F, M}”* is 
2” - 1. Let a: G R" 

n 

( 2 . 1 ) 

r—1 


where G G GG to G Z" and j G Nq. Then ■ j G Nq, G G G-’, to G Z"} is an 

orthomormal basis in L^(R”). 

Define 

■= 250^(2-^a: - k) := 2i'ip^{2^x - k) 

where fc G Z. Then, we define 

E + ^) = 2^'^ E + ^) - ^) (2-2) 


:=E^Pfe(2^ + ^) = 2^" + (2.3) 

fez fez 

Obviously, tpj^k, <j>j,k are 1-periodic functions belongs to L^([0,1]). 
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Define 

n 

^Gmi^) = 2^^ n - rUr), (2.4) 

r—1 

By Proposition 19.11 we know that ip^, tp^ included in the functions with radial 
decreasing L^-majorants, that is 

\Pj^{x)\<R^{\x\) \i;^{x)\<R2{\x\), 

where i?i and R 2 are bounded decreasing functions belongs to L^([0,oo)). Now, 
we can use Theorem 5.9 in [15] to find that : j € NU {0}, G G , m € Mj} 

with Mj = {m : m = 0,1, 2, ■ • • , 2-’ — 1} is an orthomormal basis in T^(T"). At this 
point, considering Corollary 5 in m and Definition l9.2l in the Appendix, we easily 
obtain the following theorem for wavelet characterization of variable index Besov 
and Triebel-Lizorkin space on periodic domain. 

Theorem 2.8. Let s(-) G fl and p{-) G The symbol A 

stands for B or F and so does a symbolize b or f, respectively. 

(i) Let 0 < q < 00 (p'^ < 00 in the F-case) and 

k > max(tTp — s~ ,s^) {crp,q in the F-ease), 
where = n - l) and (Jp,, = n - l) • Then f G 5'(T") 

sf 0 

belongs to ^ if, and only if, it can he represented as 

00 

f = HY^ ^ KW.F (2-5) 

1=0 GeGJ mGMj 

with Mj = {m : TO = 0,1, 2, • • • ,2-^ — 1} and the series expansion it if. 51) is uncondi¬ 
tional convergence in and in any space j ^(T"), where cr{x) < s{x) with 

inf(s(x) — cr(x)) > 0 and a{x)/s{x) —i 0 for \x\ —>■ oo. The representation i ti).,5)) is 
unique, we have 

and 

is an isomorphic map from A^j j ^(T") onto j Moreover, if in addition max(p+, q) < 

oo, then {^cmljeNo.GeG^.meMj nn unconditional basis in A®| ^ ^CIT”)- 

(ii) Let q{-) G p'°ff(T") with 0 < p~ < p'^ < oo, 0 < q~ < q^ < oo and let 

k > max((Tp_q — s“, s’*") 

with ap^q = n - l) • Then f G 5'(T”) belongs to ^p(.),,(.)(T”) if, and 

only if, it can be represented as i t A 5|) with A G /p(.j g(.)C^")) with unconditional 

convergence in <S'(T”) and in X G j ^(T”). The representation i2.5\) is unique, 
we have 

A^G™ = A^G™(/) = 2^'^ </,«'^G,n> 
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and 

is an isomorphic map from ^p^.) ■^p( j 

At the end of this section, we need to introduce the following notation which 
used frequently in the sequel. 

/ OO 

E E E ( 2 . 6 ) 

j=o GeGJ meMj 

where Xq^ are defined as in Theorem 12.81 

3. Variable order Besov prior 

For the reader’s convenience, let us recall the general setting stated in [5T] for 
our purpose. Denote J to be an index set, let denote an infinite sequence 

in the Banach space X, with norm || • ||, of M-valued functions defined on a domain 
D. In the following, for simplicity, we assume D = T” to be the n-dimensional 
torus. We normalize these functions so that ||(()j|| = 1 for j £ J. We also introduce 
another element toq G A, not necessarily normalized to 1. Define the function u by 

u = mo + E Ujfij (3.1) 

j&J 

By randomizing u := {uj}j^j we create real-valued random functions on D. (The 
extension to K”-valued random functions is straightforward, but omitted for brevity.) 
We now define the deterministic sequence 7 = {'JjjjeJ and the i.i.d. random se¬ 
quence 5 and set Uj = 'Yjfj- We assume that f is centered, i.e. that 

E(^i) = 0. Formally we see that the average value of u is then mg so that this 
element of X should be thought of as the mean function. 

In the following, we take X to be the Hilbert space 

A:=L2(t") = E : T"^K : f \u{x)\^dx<oo 

L JT" 

of real valued periodic functions in dimension n < 3 with inner product and norm 
denoted by < •, • > and || - || respectively. We then set mo = 0 and let 

J = {j = 0,1, • • • , M,} 

appeared in the Theorem 12.81 which is an orthonormal basis for X. Consequently, 
for any u G X, we have 

00 

“(*) = E E E ^GmEm with U^Gm =< Em >> (3.2) 

i=0 GeGi meMj 

where is the wavelet basis stated in the Theorem 12.81 Given a function 

w : T” —>• M and the (EmI defined in (13.21) we define the Banach space by 

E5:; = (m:T"^R : 


< OO 


(3.3) 
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with 


p*C-) 

-!(■).l(-) 


<!(■) 


l/9(-) 


(3.4) 


Lo(')(T") 

where ipjix) = (p{2~^x) and (/?(•) is a smooth decompositions of unity [HI [19]. 
By Proposition 5.4 in [T7]j we can define the space -Sp(.j p(.) appropriately which 
is equivalent to ^p(.) p(.) if s S L°°. So here, we do not need to develop a full 
theory for the space -Spp) p(.p and just understand it as ^pp) p(.)- Hence, our space 
dehned in (13.3p just the usual variable index Besov space with p(-) = g(-). 
(Although the space defined in [l^ is in the whole space R", it can be adapted to 
periodic case T”.) 

Now, as in the general setting, we assume that = Icm^Gm. where ^ = 
{^Gm}i=i, 2 .--.oo,GeGymeMj is an i.i.d. sequence and 7 = { 7 Gm}i=i, 2 ,--.oo.GeGymGM,- 
is deterministic. Here we assume that is draw from the centred measure on K 
with density proportional to exp (—i fj„ K{dy)) for some I < q~ < < 00 , 

where /«(■) is a probability measure. We refer to the measure with above density as 
a generalized (jf(-)-exponential distribution. Noting that if q is constant, it is just a 
(j-exponential distribution Hence, our generalized g(-)-exponential distribution 
is a natural extension of g-exponential distribution and include Gaussian, Laplace 
distribution as special case. For s(x) > s“ > 0 and d > 0 we define 

l/q+ 


iGr, 


_ 2 —j{s {2 ^ m)-\-nj2—njq'^ 


1) 


We now prove convergence of the series 


N 

j=0 GgGJ mGMj 




to the limit function 


EE 


E 


^Gm^ Gim 


^Gm^Gm^ 


^Gm = iGn 




^Gm = 7Gr, 




(3.5) 


(3.6) 


(3.7) 


j=0 GeGi meMj 

in an appropriate space. In order to understand the sequence of functions 
fully, we introduce the following function space: 


:= <{ u : T" X H 


3 A > 0, (Am) < 00 
«(■) 


where 


Uk = {‘Pku)'' 


with ipk defined as in (13.4L and 
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As mentioned in section \TJ\ if p^t( ) is a convex semimodular on Lp^ ^ 

Then ^ is a normed space with the Luxemburg norm given by 

II«IIl 2»(0;B*«,) = inf > 0 : < l| ■ (3.9) 

In order to keep the fluency of our statement, we list the proof of is a convex 

-®9() 

semimodular in section[5] Now, we clarify the relation for our space Lp^'\fl] 
with the usual constant q,t space Lp{n;B^) used in [53]. Let q, t in (13.81) to be 
constants, we have 


p OO n 

PBt{'^)= / / \uk\‘^dx¥[duj) 

" k=o 


= E(hllL )■ 

Hence, our variable space indeed is a natural generalization of the usual space 
Lp(H;i?‘). Define 


Kn ^ b - {A^^ljgNo.GeGL 


rEs(.) < OO ^ , 
g(-) 


(3.10) 


where 


i(-) 


i, i=0 GgGJ meMj 



\ l/9(-) 




L9(-)(T'*) 


and 


\ / Jjn 


With these definitions, before going to our one main results in this section, we need 
the following Lemma, which is proved in the Section |8l 


Lemma 3.1. Let s(-) S n (7;°^(T”) and q{-) e 7’'°9(T"). Let 

k > max((Tq — s~, s'*'), 

where aq = n . Then f G <S'(T”) belongs to Lp^ ^(D; -Sq[.]) if, OL'nd 

only if, it can he represented as 

OO 

f = HY^ H WlthXGb^^f\ (3.11) 

i=0 GeGi meMj 

with Mj = {m : m = 0,1, 2, • • • , 2-^ — 1} and the series expansion is uncon¬ 
ditional convergence inS'iMT). The representation is unique, we have 

A^G™ = A^G™(/) = 2^'^ </,«'^G,n> 

and 




is an isomorphic map from Lp^ ^(D; onto 6^^) 




bEs{-) 
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Now we can prove the following theorem, which gives a sufficient condition, on 
t(-), for existence of the limiting random function. 

Theorem 3.2. For t,s€ n L°°(T"), q S and 

sup ( t{x) — s(x) + -^ ) < 0 

xeT" V 

the sequence of functions given by ifg.dl) and I13.5\) with draw from 

a centered generalized q{-)-exponential distribution, is Cauchy in the Banach space 

Lp^ ^ ■ Thus the infinite series jS. 7| ) exists as an a limit in space Lp^ ^ 

for all sup^.gT" (tix) - s{x) + < 0. 

Pvoof. By (j3.7p cLiici tli6 LcninicL ixn we obtain 

For M > N, every A > 0, we have the following estimate 

M 


(3.12) 




=N+l GgG^ meMj 

» M 

<C'max(A«^,A«')(5-i / E E E 


j=Ar+l GgGJ mGMj 


M 


<C'max(A« ,A« )(5"1 Y 2^"^ ^m)-s {2 ^■m.)+n/q^)^-jr. 

j—N-\-l mGMj 


(3.13) 


where we used 


E 




< c 


< c 


/Rn{| 5 |>i} 
+ C ; 

/Rn{|5|>i} 

+ C 


ier%xp (^P-j^j^\‘ii-)n[dx)di 


exp ( -- / |C|«(")«(dx)de 


/Rn{|c|<i} 

|^|«%xp(-^|^r'd^ 


/T" 


1^1’ exp () < oo. 


./Rn{|4|<i} 

The sum on the last line of (13.131) tends to 0 as iV —)► oo, provided 

sup ( t{x) — s(a;) + E ) <0. 

a:GT" V 9^ / 

Hence, by Lemma 2.1.9. in [20], we obtain 


^lini^||{2^t(u^^-u^i)}||^..(.) 


= 0 . 


Finally, by (I3.12L we complete the proof. 


□ 
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Remark 3.3. Here we give an intuitive meaning for the random series which we 
defined in (1X71) . Assume the probability measure /c(-) in the centered generalized 
( 7 (-)-exponential distribution is an uniform measure that is K{dx) = dx in T". Since 
orthonormal basis and 

OO 

“ = Z! Z (3.14) 

j =0 GeG^ meMj 


with = 2 ^m)+n /2 n/q+) denote we have 

7 / 


n n n 

i=0 GeG^ meMj 

OO 

=n n n 

j^o GgGJ mGMj 

OO 

=n n n ^^p 

j^o GgG 3 mGMj 

OO 

n n e^p 

i=0 GeGJ meMj 

1 


l7i 


J \-Qi^)\nd 


Gm I 


^Gm 




<exp --min{(57F,5}p .(.)(m) 

V 2 ^q(,) 

Thus, informally, the Lebesgue density of u can be controlled by a Lebesgue density 
proportional to exp i —h min{(57^, S}p„i,(.) (u) ). Since p„s(.) (u) related to the space 


Jjn 


fTheorem I2.8I1 . and the space is a generalization of constant 

index space, we may guess the Lebesgue density of u related to a Lebesgue density 

similar to exp ( —-1 minl^T^^, (5}p„s( ) (u) ). At least, if we choose q as constant 
. A . . ’(■) / 

and just s is a function, we will have an equality which informally means the 
Lebesgue density of u is proportional to exp iu)j ■ Hence, the probability 

measure we defined may be related to the space i?*| |(T"). So we may say that u 
is distributed according to an measure with parameter 5, or, briefly, a 

((5, measure. 


K-) 


g(-) 


Theorem 3.4. Assume that u is given by Jg. 7[ ) and i3.5\) with for every 
{j, G, m} draw from a centered generalized q{-)-exponential distribution with n^dx) = 
dx that is to say «;(•) is an uniform probability measure. In other words, u be dis¬ 
tributed according to a (6, .Bg(.j) measure. In addition, we assume t,s G (7*°®(T") D 
L°°(T”), q G )P*°s(T”) and t{x) — s{x) + ^ 7 ^ 0 for every x G T”. Then the 
following are equivalent: 


(1) Pgt(){u) < OO P-a.s.,- 


q() 


(2) E ( exp ( ap^tc ) (u) 


< OO for any a G [0,(5/2); 
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(3) sup^gT" {t{x) - s{x) + < 0. 

Proof. (3) ^ (2). 

Since sup^g^” (t{x) — s(a;) + < 0, there exists a negative constant /3 < 0 

such that sup^gjn (t{x) - s{x) + ^) < /3 < 0. Let if to be a large enough positive 
constant, then we have 

E I exp ( apgto) (u) 


<?(■) 


= n n E(exp(ari f 

i=0 GeGi \ mGM,- , 


<n n E|exp|ari ^ /_ 

j=0 GeGi 


n 

j=0 GeGi 


4 exp (aJ-i 4„ - i 4„ 




< 


n 

u=0 


4 exp {a6-^ 4„ ^ rfg 

/Rexp(-HT" l^l'^'^^^dx) d^ 


K 


If we want to prove the above infinite product converge, we only need to prove the 
following summation converge |26| . 


f2 J (^exp J - 1^ exp J d^ 


4.g g ^ il »P (4 X. <■ 

Now we concentrate on the first summation term in the above integral. 

oo oo , 


(3.15) 


EE 

j—0 k—l 


{ad 


kl 






fc = l 


/ 


l/k 


fc=i • \ A" \j=o I 




< 


< 


E 


(a(5 


1 - 2^*?" ^ k\ 

fc=i 




^-^exp 
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Substituting the above inequality into (I3.15L we obtain that 

J ^exp y — 1^ exp i J 

where we used a5~^ < 

(2) => (!)• 

If (1) does not hold, p„t( ) (u) is positive infinite on a set of positive measure S. 
Then, since for a > 0, exp I ap^t( ) {u) J = +00 if p^to) {u) = + 00 , and 

V ^<z() / '°«() 


E ^exp ^ap^to) (m)^ ^ > E ^ s exp ^ap^to) (u) 


we get a contradiction. 

(1) => (3). 

Since PoU ) (u) < 00 , we easily know 





j=0 GeGi meMj 


Hence, for almost all x G T", the integrand in the above formula is finite. Choose 
X G T" such that 


j=0 G€G^ meMj 


So for every j, there is m = mxj such that 

00 

j=0 GgGJ 

If t{2~^mx,j) — s{2~^mx,j) + n/q^ > 0, we have 


E E 1^' 

i=o GeG3 


Grrix 


\l{^) 


< 00 . 


Since there exist c,C such that 0 < c < < C < 00 for every 

{j,G,m}, this contradicts the law of large numbers e.g. Theorem 19.31 So we 
obtain t{2~^mxj) — s(2~^mx^j) + n/q'^ < 0 for infinite j. By the definition of 
we know 2~^mxj —>■ x, hence, we find that 

Tl 

t(x) — s(x) H—V < 0. 

In addition, by our assumption and the continuity of t(-) and s(-), we finally obtain 

Tl 

t[x) — s(x) H—ij: < 0 

for every x G T”. □ 
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Remark 3.5. Theorem 13.41 assumed t{x) — s(x) + ^ 7 ^ 0 and k{-) is a uniform 
probability distribution which seems technical, however, for the constant q, t, s case 
these conditions are all satisfied naturally. How to remove these conditions will be 
left to our future work. 


In the previous two theorems, we proved basic properties for random variables 
construct from infinite series (EH). Now we study the situation where the family 
have a uniform Holder exponent a and study the implications for Holder 
continuity of the random function u. We assume that there are C,a,b > 0 and 
a G (0,1] such that, for all j > 0, 

xGT", 

We also assume that a > 6 as in Ha- 

Theorem 3.6. Assume that u is given by {3. 7[ 1 and iS.5\} with draw from a 
centered generalized q{-)-exponential distribution. Suppose also that h3.16\) hold and 

that s G (T") n q G inf^eT- s{x) > n (b + ^ + \e{a - 5)) 

for some 9 G (0, 2). Then P-a.s. we have u G for all ft < 

Proof. Here we need to use Theorem 19.41 listed in the Appendix which is a variant 
of the Kolmogorov continuity theorem. Denoted as in the Theorem l9.4l but use our 
series (ITtI) . we obtain 

OO 

= E E E 

j=0 GeG3 meMj 

i=0 GeGJ meMj 


< C 


OO . . . OO 


3=0 


3=0 


and 


S2<C'^2 E E iT'Cml^ ^W^GrriWl^^llGml 


2^jna6 


where 


j=0 GeG^ meMj 

OO , V 

< = c'^2 

3=0 3=0 


2s(2-^'") 2 „ 

Cl = —- - -- - 25 > 0, 


-3nc2 


and 


C2 = 


2s(2-J”") 


- — - 2b - 9{a - b) > 0. 


n q~' 

We need ci > 0 and C 2 > 0 and by our assumption a > 5, we only need C 2 > 0. 
Our assumption infa;gT" s{x) > n (^b + -^ + ^9{a — b)^ just assure C 2 > 0. So by 
Theorem 19.41 we can conclude our proof. □ 
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Remark 3.7. If let the mean function is nonzero and satisfies 


|wo(a;)| < C, X £ D, 

|mo(a;) - mo{y)\ < C\x - y|“, x,y £ D. 

Then the result of Theorem 13.61 still holds. 

Theorem 3.8. Assume that u is given by ({S. 9| ) and \3.5\) with drawn from 
a centered generalized q{-)-exponential distribution. Suppose also that 'I'g^, with 
{j = 0,1,-- - ,oo,G £ G^m £ Mj} form a basis for with t~ > 0, q{-) £ 
and t £ (T") n L°“(T"). Then for any 

sup ( t{x) - s(x) + -^ ) <0, 
xejn V 9 / 

we have u £ C'*^'^(T") V-almost surely. 


Proof. For any k > 1, using the definition of PnH ) (u), "we can write 

9 () 


P f>^() (l^) — Gs^rn 


i=0 GeGJ meMj 


For every k £N there exists constants Cm, Gm such that 0 < Cm < E — 

Gm < oo. Since each term of the above series is measurable we can swap the sum 
and the integration and obtain 


E 


V ''Of) ) 


< oo. 


From the above inequality, we obtain that p„t( ) (u) < oo P-almost surely. So we 

) 

know that ||M||^t(.) < oo P-almost surely. Since 


sup ( t(x) - s(x) + E I <0, 
a:GT" V 9+/ 

we can choose fc large enough SO that < s{x) —-^^—t{x). Then the embedding 

^ C'‘(')(T") [H] for any G satisfying t{x) + < G{x) < s{x) - ^ 

implies that ||it||( 7 t( )(T") < cxd P-almost surely. It follows that u £ C*^'^(T”) P- 
almost surely. □ 


Remark 3.9. If the mean function mo belongs to the result of the above 

theorem holds for a random series with nonzero mean function as well. 


4. Bayesian approach to inverse problems 

In this section, we generalize results in m and [23] to our setting. Let X, Y 
be separable Banach spaces, equipped with the Borel cr-algebra, and Q : X ^ Y 
a measurable mapping. We wish to solve the inverse problem of finding u from y 
where 


y = G{u) + p 


(4.1) 
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and T] G Y denotes noise. Employing Bayesian approach to this problem, we let 
(u, y) G X X Y he a. random variable and compute u\y. We specify the random 
variable (u, y) as follows: 

• Prior: m ~ yo measure on X. 

• Noise: y ~ Qo measure on Y, and y -L u. 

The random variable y|u is then distributed according to the measure the 
translate of Qo by Q{u). We assume throughout the following that Qu ^ Qo for u 
fiQ-a.s. Thus we may define some potential $ : X x T R so that 

^(u) =exp(-$(u;y)), (4.2) 

and 

J exp(-$(u;y))Qo(c^ 2 /) = 1- (4.3) 

In the previous section, we construct probability measure yo which are supported 
on a given variable order Besov space Following, we show how to use of such 

priors yo may be combined with properties of <i>, defined above, to deduce the 
existence of a well-posed Bayesian inverse problem. We firstly list the following 
conditions of <i>: 

Assumptions 1: Let X and Y be Banach spaces. The function $ : X x F —R 
satisfies: 

(i) there is an oi > 0 and for every r > 0, an M G R, such that for all u G X, 
and for all y G F such that ||y||v < r, 

‘h('u;y) > M - ai||u||x; 

(ii) for every r > 0 there exists K = K{r) > 0 such that for all u G X, y G F 
with max{||tt||xj \\y\\Y}<r 

$(u;y) < K; 

(hi) for every r > 0 there exists L = L{r) > 0 such that for all Ui,U 2 G X and 
y G F with max{||ui||x, ||m 2 ||a:, ||2/||v} < r 

|$(Mi;y) - $(u 2 ;y)| < T||ui - uzHx; 

(iv) there is an 0:2 > 0 and for every r > 0 a C G R such that for all yi, y 2 G F 
with max{||yi||v, ||y 2 ||v} < r and for every u G X 

l^'(M,yi) - $(u,y2)| < exp(a2||'«||x TC') Wvi -y2||y- 

Now we can give the following theorem for well-defined problem. 

Theorem 4.1. Let <i> satisfy 14-^ Assumptions 1 (i)-(iii). Suppose that for 
some q G P^°®(T”), t G nL°“(T”), is continuously embedded in X. 

There exists <5* > 0 such that if po is a (d, j) measure with 

sup ( t{x) — s{x) + -^ ) < 0 

xeT" V r J 

and 5 > 5*, then is absolutely eontinuous with respect to po and satisfies 

dp^ ( ^ 1 

dpo Z{y) 


exp(-$(u;y)) 


(4.4) 
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with the normalizing factor 

Z(y)= exp(-$(M;y))/7o(du) < oo. (4.5) 

Jx 

The constants* = 2max{c® where Ce is the embedding constant satisfying 

||m|U < CeWuWgtO- 


Proof. Define 7ro(du, dy) = fiQ{du) ® Qo(dy) and '!T{du,dy) = yLo{du)Qu{dy). As¬ 
sumption 1 (iii) gives continuity of $ on A and since /io(A) = 1 we have that 
$ : A —i> R is /iQ-measurable. Therefore tt ttq and tt has Rando-Nikodym de¬ 
rivative given by (14.21) . This then by Theorem 6.29 of [5S], implies that yT{du) is 
absolutely continuous with respect to yo{du). This same lemma also gives (14.41) 
provided that the normalization constant (14.51) is positive, which we now establish. 
Since = 1, we note that all the integrals over A may be replaced by inte¬ 
grals over for any sup 3 ,gT" (^( 2 ;) — s{x) + < 0. By Assumption 1 (i), we 

note that there is M = M (y) such that 



exp {-^{u;y)) fj.o{du) 


exp(ai||it||x - M) yo{du) 


By Lemma 3.2.5 of EQ], we know that 


when Pgt( ) (u) > 1 or 


1 ^ <ct 


PnHoiu) 


<j(-) 


\x <4 


<,(■) 


when p„t(.) (it) < 1. Hence, if ||it||A > 1 then we have 


«(■) 


Z{,y)= / exp(-$(ii;y))/io(dii) 

g(-) 


< 


g(-) 


exp ( «! max{c® , c® }y„t( ) (u) — M \ yo(dit). 


■j(-) 


This upper bound is finite by Theorem 13.41 since S > 2ai maxjc^ ,c|^}. For 
llitllx < 1, we have 



exp(-$(u;y))/io(dit) 


< 



exp(Qfi — M)/io(dit) < 00 . 


Now, we prove that the normalisation constant does not vanish. Let R = K(p„t( ) (it)), 

9( ) 

we know that R G (0, 00 ). As p^to (it) is a nonegative random variable, we obtain 
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that ^o(pr^t( ) (u) < R) > 0. Taking r = max{ y Assumption 1 (ii) ensure 

9() 

that 


z{y) 



exp{-^{u\y))po{du) 


> / exp{—K)fio{du) 

= exp(-Ar)yo(PRt() (u) < R) 


which is positive. 


□ 


We now show the well-posedness of the posterior measure with respect to the 
data y. Recall that the Hellinger metric dneii is defined by 


dHellifJ-, d-') 





dv, 


where v is reference measure and with respect to which both p and p' are absolutely 
continuous. The Hellinger metric is independent of the choice of the reference 
measure v. For a review of probability metrics we refer to |29) . The following 
theorem can be proved by using similar ideas from Theorem 3.3 of |23| . The minor 
differences are that we need to use p„t( ) M instead of \ when we need to bound 

||m||x- The same situation appeared in the proof of Theorem 14.11 so we here omit 
the details of the proof for concisely. 


Theorem 4.2. Let <i> satisfy i4-d\ l and Assumptions 1 (i)-(iv). Suppose that for 
some q G t G C'jof (T”) fl L°°(T”), is continuously embedded in X. 

There exists <5* > 0 such that if po is a ((5, j) measure with 

sup ( t{x) - s{x) + ^ ) < 0 

V r J 

and S > 6*, then 

dHeu{p^,p^ ) < C\\y - y'Wv 

whereC = C{r) withmacx{\\y\\Y, Hy^Hv} < r. The constants* = 2max{c® ,c|^}(q;i + 
202 ), where Ce is the embedding constant satisfying ||u||x < Ce||u|| „t( ) • 


In the last part of this section, we consider the approximation of the posterior. 
Consider to be an approximation of $. Here we state a result which quantifies 
the effect of this approximation in the posterior measure in terms of the approxi¬ 
mation error in 4>. 


Define p"^'^ by 


and 


dpy 


N 


dpo 


-{u) = 


Z^{y) 


exp , 


rN 


(y) = f exp (-4>^(u)) po{du). 
Jx 


(4.6) 


(4.7) 


We suppress the dependence of $ and 4)^ on y here as it is considered fixed. 
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Theorem 4.3. Assume that the measures ^ and are both absolutely continuous 
with respect to fio, and given by o and respectively. Suppose that <I> and 

4*^ satisfy Assumption 1 (i) and (ii), uniformly in N, and that there exist 03 > 0 
and C S K such that 

|4>(m) - < exp(a3||M||x + C)(p{N), 

where (p{N) —)► 0 as TV —>• 00 . Suppose that for some t S H L°°(T”), 

q g is continuously embedded in X. Let ^0 be a measure 

with 

sup ( t{x) — s{x) + -^ ) < 0 
xeT" V 9'^/ 

and 6 > 2 max{c| ,c^^}(q;i + 203 ) where Cg is the embedding constant satisfying 
||m||a: < Ce||M||^t(.). Then there exists a constant independent of N such that 

dHell{pi,li^)<Cg,{N). 

The proof of the above theorem similar to the proof of Theorem 3.3 of [55], 
the differences and difficulties can be overcome by same idea used in the proof of 
Theorem 14.11 


5. Variational methods 

MAP estimator in the Bayesian statistics literature ISO] is an important concept, 
it specifies the relationship between Bayesian approach and classical regularization 
technique. It is well known that for non-gaussian prior we can hardly obtain a 
rigorous relation between prior measure and regularization term in inhnite dimen¬ 
sions. Even for a simpler constant index Besov prior there is no complete theory 
[23] . Only recently, Tapio Helin and Martin Burger m addressed this issue in 
some sense. Here, as stated in Remark 13.31 we can get an upper bound for the 
probability density, and for constant q case we can get an equality. So we faced 
more complex situation compared with constant Besov case and in this section we 
just give some partial illustration. 

We define the following functional 

I{u) = $(m) -b ) (u). (5.1) 

Z q(-) 

Intuitively, the minimizers of the above functional or some variant of (ED) may 
have highest probability measure when a small balls centred on such minimizers. 
For more explanations, we refer to the Gaussian case [32]. For this functional we 
have 

Theorem 5.1. Let Assumption 1 (i),(ii) hold, assume that poiX) = 1. Let s S 
nL°“(T”), q G I < q~ < q{x) < < 00 and compactly 

embedded in X . Then there exists u G -®q(.) such that 

I{u) = / := inf |/(u) : u G | . 

Furthermore, if {un} is a minimizing sequence satisfying I{un) —S' I{u) then there 
is a subsequence {un'} that converges strongly to u in B^^^ y 
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Before proving this theorem, we need to prove the following lemma which is of 
independent interests. 

Lemma 5.2. Let s S (7'°^(T") nL°°(T"), q e P'°9(T"), 1 < g” < q{x) <q+ <oo, 
then the dual space of is where 

— + ^— - 1 
q{x) q{x)’ 


Proof. Step 1. Let s G Cl°^nL°°, q G We prove in this step that (R") C 

where stands for the dual space of Let / G 

, define fk = J2l=-i{Pk+rfY where (pi defined as in (13.41) with is a 
smooth decompositions of unity [niiis]- Then we know that 

OO 

f = Y^{pkfkr iii5'(R") (5.2) 

fc=0 


and 


l|2 Vfc|lL9'( )(R",f<l'(-)) 
= inf < A > 0 : 


/R" 


e^—js{x)q' (x) 

— JTiP) — \iPjfk:yix)\‘>'^^^dx < 1 


j=o 


<Cmii ^ >0: 


<C\\f\\s-H-G 

<l(-) 


i^—js{x)q' (x) 

3=0 (. 2 ) 


(5.3) 


Take <p G 5(]R"), then we have 


\fip)\ = 


00 1 

EE (V^fc+r/fc+r) )) 


k—0r—-l 
00 1 


± p 

= ^ ^ / 2-<-'>^Mx)r^^^^p{ipm+rP'')dx 

< C||2 ^^/fc|lL<!'(-)(Rn,f<i'( ))l|2®^ '^^iF{(pkPk+rP'^)\\L<i(-)(S.^/'i( )) 


Hence, now we proved i?~^®|’^(IR”) C i?®|’j(]R”)'. 

Step 2. In this step, we need to prove i?*|’j(R")' C (R"). Since / G 

i3®|'j(R") —)• i® 0 . one-to-one mapping from i3^|'j(R") onto a 

subspace of L«(')(R", f'J^)), every functional g G (B^^'j)' can be interpreted as a 
functional on that subspace. By the Hahn-Banach theorem, g can be extended to 
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a continuous linear functional on where the norm of g is preserved. 

If G 5(K”), considering Corollary 1 of Theorem 8 in Chapter 13 of [42] we have 

p OO 

9i‘P)= '^gk{x)[ipk‘fY {,x)dx (5.4) 

fc =0 

where ^<,'(.)) is equivalent to the operator norm of g. (|5.4D implies 

that 

OO 

g = ^T{tpkgl) in 5'. (5.5) 

r\nk 

Define r]km{x) = (i_|_ 2 fc| 2 :|)m with m is a large enough constant, since ip can be 
chosen to be a radial smooth function with compact support, we know that it can 
be controlled as follows 


2^^p^{2^x)<Cgkm{x). 


(5.6) 


Now, we have 




inf|A>0: f 1 

I JR" A J 

Jr" A9 J 


= inf A > 0 

<(7112 '^’^Pk * '?fc|lL<?'(-)(R",^<i'( )) 

<(7||2 * 5fc|lL<!'(-)(R",.f9'{-)) 

<(7|12 **■ ^^5fc||i9'(.)(Rn_^,'(-)) 

where the last inequality follows from Lemma 5.4 of [24]. (Here we need a small 
modihcation of Lemma 5.4 of [24] . however, the modification is straightforward so 
we omit for concise.) With the above estimates, obviously, the proof is completed. 

□ 


Although the above proof is for the whole space R”, it is also valid for periodic 
domain T". Now, let us come back to the proof of Theorem 15.II 


Proof. For any <5 > 0 there is A^ = A^i(i5) such that 

I < I{un) <1 + 5, \/n> Ni. 


Thus 


:rPB«()(Mn) < M- 
Z g(-) 


\/n> Ni. 


The sequence {un} is bounded in By the above Lemma [521 we know that 

j is reflexive, there exists u € such that m„ ^ it in B^^^ y By the compact 
embedding of B^^ !^ in X we deduce that —>■ u, strongly in X. Noticing that 

) (u) is lower semi-continuous by Theorem 2.2.8 of [20]. Now, we can use similar 

g(-) 

ideas of Theorem 2.7 in m to complete the proof. □ 
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6. Application to backward diffusion problem 

In this section, we apply our theory to inverse source problems for integer order 
diffusion equation and fractional order diffusion equation. 


6.1. Integer order diffusion equation. For simplicity we consider periodic do¬ 
main T". Define the operator A as follows: 

A = -A, V(A) = ff^(r). 

Consider the diffusion equation on T" with periodic boundary condition, writing it 
as an ordinary differential equation in JJ: 

+ Av = 0, v(0) = u. (6.1) 

Define G(u) = e~'^u, ^ to be an operator defined as follows 

i{G{u)) = {G{u){xi), G{u){x2), ■■■ , G{u){xk))^ ( 6 . 2 ) 

where AT is a fixed constant. Then we have the relationship 

y = e{G{u))+ri (6.3) 

where rj = is a mean zero Gaussian with covariance T, y = is the 

data which we measured. Now we can show the well-definedness of the posterior 
measure and its continuity with respect to the data for the above inverse diffusion 
problem. 


Theorem 6.1. Consider the inverse problem for finding u from noisy observations 
of G{u) = u(l, •) in the form of S6.S\} . Let y,o to be distributed as a variable index 
Besov prior (d, with s £ (^^"^(T") nL°°(T"), q G 'P*°®(T”), infa;gT" s(®) > ^ 

and (5 > 4. Then the measure yT{du) is absolutely continuous with respect to /ig 
with Radon-Nikodym derivative satisfying 


dpy 

-i—w 

dfio 


^^exp(-$(M;2/)) 


where 

^{u-y) = i|r-V2(y _£(g(u)))| 2 _ l|r-V2y|2 


and 

^ X*'-’ “ ^(^( m )))^ + Mdu) 

g(-) 

with sup,j,gT" — s{x) + < 0. Furthermore, the posterior measure is con¬ 

tinuous in the Hellinger metric with respect to the data 

dHeii{tJy,y^ ) < G\y-y'\. 


Proof. Firstly, we prove two properties about the operator £{G{-)). 



BAYESIAN INVERSE PROBLEMS 


23 


Property 1: For large enough £ > 0 and a small constant e > 0, by Sobolev 
embedding theorem, we have 

|£(G(n))| <iF||G(n)|U=o = ||e-^n|U=o 

< \\A^e~^u\\L2 + ||e“^M||i2 

<G||w|| (6-4) 

^2,2 

< G||M||^t(.) 

<l(-) 

where we used the fact that ^ A > 0, is a bounded linear operator from 2 

to i ?2 2 i ® we also used embedding theorems for variable index 

Besov space m- 

Property 2: Let u\, U 2 be two different initial data for heat equations. By 
similar idea from the proof of Property 1, we have 

|£(G(ni)) - £(G(u2))| < K\\G[ui) - G(ir2)IU» 

= ^l|G(ui-M2 )||l=o (0 5) 

< G||mi -M2||^t(.). 

g(-) 

Now let X = £3*1 j. With the above Property 1 and Property 2, it is straightforward 
to know that satisfy Assumption 1 (i)-(iv) with oi = 0 and 02 = 1- By 

Theorem 14.11 and Theorem 14.21 we immediately obtain our desired results. □ 


6.2. Fractional order diffusion equation. For fractional order diffusion equa¬ 
tions, there are numerous literature. For the well-posedness theory, we refer to 
[33l [Ml ISa [36]. Here we just treat fractional diffusion equations on periodic do¬ 
main as follows 

dfv{t, x) + (-Afv(t, x) = 0, t>0, X e T", 

v(0) = u, 

where 0 < a < 1 and 0 < /3 < 1 and 9“ stands for Caputo derivarive of a order, it 
can be defined as follows 

'■= 1 [ (.t- s)~°‘fis)ds, 

r(l - a) Jq 

where F(-) is the usual Gamma function. Define the operator A as follows: 

iL=(L2(T"),<.,.>,||.||) 

A = (-A)^, P(A) = ij2/3(T"). 

Consider the heat conduction equation on T" with periodic boundary condition, 
writing it as an ordinary differential equation in H: 

dfV + Av = 0, v{0) = u. (6.7) 

If A is a bounded operator e.g. a positive number, then the solution of the above 
equation (16.71) has the following form 

v{t) = Ea{—At°‘)u, 
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where Ea{-) is Mittag-Leffer function defined as 


Eciz) = ^ 

k=0 


r{ak + 1) 


For more properties about Mittag-Leffer function we refer to [37l |38] . In [39] , they 
proposed fractional operator semigroup which characterize the solution of abstract 
fractional cauchy problem (IQ) . Since operator A in our setting can generate a 
fractional operator semigroup, we can define G{u) = Ea{—A)u, i to be an operator 
defined as follows 


t{G{u)) = (G(m)(xi), G{u){xi), • • • , G{v)[xKyf (6.8) 


where K is a fixed constant. Then we have the relationship 


y = i{G{u))+r] 


(6.9) 


where rj = is a mean zero Gaussian with covariance T, y = is the 

data which we measured. Reviewing the proof of Theorem 16.11 the key points are 
the estimates about the operator £(G(-)). Here, we meet more crucial situations 
for the fractional diffusion equations have no strong smoothing effect as normal 
diffusion equations. For a more complete illustration, we refer to [13]. Mittag- 
Leffer function has only polynomial decay rate which restrict the smoothing effect. 
More precisely, we list the following decay rate estimates 


Lemma 6.2. 


If Q < a < 2, y, is an arbitrary real number such that 


na 


< y < min{7r,7ra}, 

then for an arbitrary integer p > 1, when \z\ —>■ oo the following expansion holds: 


^-k 


a 1 (p — 


0(|z 




where \arg{z)\ < y. 


Based on the above observation, we must restrict the fractional order a, /3 to 
some appropriate interval to gain enough smoothing effect to obtain the forward 
operator is lipschitz continuous. More precisely, we can obtain the following result. 


Theorem 6.3. Consider the inverse problem for finding u from noisy observations 
of G{u) = u(l,-) in the form of 116.9\} with 0 < a < 1 and f < /? < 1. Let 
X = L^(T"), yo to be distributed as a variable index Besov prior (i5, j) with 
s € G;°®(T”) n L°“(T"), q G P^°9(T"), infj;^!" s(a;) > ^ o,nd 5 > A. Assume 
t G G;°f (T”) n L°“(T") and 


n 

q~ 


n , 

- ^ <t~ <t+ <s- 


n 



Then the measure y^{du) is absolutely continuous with respect to yo with Radon- 
Nikodym derivative satisfying 


dyy 


1 
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where 

^{u-y) = i|r-V2(y _£(g(u)))| 2 _ i|r-V2y|2 

and 

Z{y) = Jx exp - ^(G(u)))p + ^iQ{du). 

Furthermore, the posterior measure is continuous in the Hellinger metric with re¬ 
spect to the data 

dHeii{yy, ) < C\y-y'\. 

Proof. In order to prove the above theorem, we first give the following estimates. 
Denote fe{i € Z’") to be the fourier coefficient of function /, then we have 

^ (|£|2/3£;^(_|£|2/3))2|j^|2 

iez^ ( 6 . 10 ) 

^ C'||/|||2(Tr.), 

where we use Lemma 16.21 Considering ^ — § < t~ < < s“ — for an 

arbitrary small positive number e > 0, we know that 

K\] ^ ^ X. 

Using (I6.10L we easily have 

|£(G(u))| < K\\G{u)\\lo^ = ||U„(-Al)u|U=o 

< \\AE^{-A)u\\l^ + \\E^{-A)u\\l^ (6.11) 

< C'I|u||l2, 

where we used ,5 > ^ to obtain the first inequality. Similarly, we can obtain 
|£(G(ui)) - ^(G(u 2))| < G||ui - U 2 IU 2 . 

At this stage, we can complete the proof as in the integer case easily. □ 

7. Conclusion 

In this paper, we firstly use the wavelet representations for function space on 
periodic domain to construct a probability measure called (<5, measure. It can 
roughly be seen as a counterpart of variable regularization terms. Using the new 
non-gaussian measure as our priori measure, we establish “well-posendess” theory 
for inverse problem as did in m- Through our study, we give another choice for 
the priori measure except Gaussian and Besov priori measure and in addition, it 
may give an another understanding of variable order space regularization terms. 

Secondly, we use our theory to integer order backward diffusion problems and 
fractional order backward diffusion problems. Especially, we prove that for time 
derivative 0 < a < 1 and space derivative ^ < /3 < 2 (n is space dimension) , the 
fractional order backward diffusion problem is “well-posedness” under Bayesian 
inverse framework. This study also reflects that the fractional order problems is 
not a straight generalization of integer order problems. When we consider fractional 
order problems, we must notice that sometimes the fractional order equations have 
totally different regularization properties compare to integer order equations. 
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8. Technical lemmas 

8.1. Properties of P^t() appeared in Section [3l 


■(■) 

<!(■) 


Lemma 8.1. Let 1 < q~ < q{-) < q~^ < oo and t{-) G C(T”). Then p^te) a 

^i(.) 

modular and continuous. 

Proof. Properties (1) and (2) in Definition 12.II are obviously satisfied. To prove (3), 
we suppose that 

Pb* . , (-^w) = 0 
«(■) 

for all A > 0. Clearly, for some ko 

[ [ {Xu) = 0 

Jn 7t" 

Since I < q~ < q{-) < q~^ < oo and t{-) G (^(T”), we easily obtain that Ukg = 0. 
Hence, we obtain it = 0. Let ^ 1, we need to prove p^to (Am) —!• p^to (m). Fix 


B ,, 
q(-) 

e > 0, choose > 0, let /i < 1 and close to 1 enough such that 

N 


<!(■) 


Jn JT" n 


< Pr 4( ) + 2e. 


B ,, 

Hence, we find p^t, ^ (Am) is left continuous with respect to A. Similar method can 

^9(0 

□ 


give the right continuous. 

Lemma 8.2. Let q €V, then p^to convex. 


Proof. Let 9 G (0,1), then 

n n OO 

P^M.)(0/+ (!-%)= / / ^2'=‘(“)«(")|0/fc(x,a;) + (l-%fc(x,w)r(")dxP(du;) 

•>(■) Jn Jt" r—1 


„E iof\ , „E 


<PBto){Qf) +Pgte){{^-^)9) 

li) <z() 

<^PRt()(/) + (l-^)PRt()(5)- 


g(-) 


iC-) 


□ 


In our case the parameter p(-) in [17] is equal to g(-), so we only need 1 < <? 
not 2 < q~ indicated by Theorem 3.6 in m 

8.2. Proof of Lemma 13.11 Concerning the proof of Lemma [3Tl we give the fol¬ 
lowing two important Lemmas. 

Lemma 8.3. Let 1 < q~ < q{-) < q'^ < oo, d > 0. For any sequence {pjffLQ of 
nonnegative measurable functions on T" denote 

Gj{x,u}) = ^2~^’^~^^^gk{x,uj). 
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Then 






where 




( OO 

j=0 


E ) 

l/9(-) 


L 9 (')(T") 


and 


Proof. It is obviously that we only need to give the following estimates 

l/q{x) , V 1/9(2:) 






< 


' i=o i=o 

\ 1 / 9 ( 2 :) 

( 00 ^ 

<C 

V/=o 


□ 


Lemma 8.4. Let q{-) £ C^°®(T”) with I < q < q{-) < q~^ < 00 . Then the 
inequality 

||{l7i;,fi * /iljeNo II ^ ll{/i}i6No llL9(-)(f^('>) 

holds for every sequence {fj{x,uj)}j^^g of Lj^^-functions for variable x and P- 
measurable functions for variable oj . 


Proof. The proof of this lemma similar to the proof Lemma 5.4 in [21], so here we 
only give the different part. Let Vi stands for all dyadic cubes with side length 2“® 
and ■ As in [24], we need the following estimate 


n OO n 
/ ^ / |77j.m * 

’•f" V=0 

p OO I p OO » 

^ E / E x3Qix)MQmduj)] dx 

JT^ „_n \ Jn .-n / 


OO « 00 




^ X 3 q(x)(Mq(|M«(")/«')) ¥{duj)dx + 

i/=0 j=o ‘ “ 

=I + II, 


II 
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where II is exactly the same as in the proof Lemma 5.4 in [24]. Next, we only give 
estimate for term I. 

/ CXD n 

E / E X3Qi^)dx 

i'=od^ j=0 QeI5„_3 

p OO 

<C 'VE dx < OO. 

With these estimates, it is easy to recover the whole proof. □ 

With the two Lemmas in hand, following the proof of Theorem 4.4, Theorem 4.5, 
we can get the local means characterizations of Lp^ ^(f2; j) by using our Lemma 
I8.3l and Lemma [8N| instead of Lemma 4.2 and Lemma 4.3 in [S^. Replacing Lemma 
5 and Lemma 9 in |16] by our Lemma 18.31 and Lemma l8.41 we can mimic the proof 
of Corollary 2 and Corollary 3 in m to give the proof of Lemma 13.11 Considering 
the proof is so long and have no new ingredient except Lemma 18.31 and Lemma l8.41 
we omit it here. 


9. Appendix 


Proposition 9.1. [16] (i) There are a real scaling function ipp S 5(]R) and a real 
associated wavelet cpM & iS(R) such that their Fourier transforms have compact 
supports, <,5 _f( 0) = 1 and 


supptpM C 


8 

2 

U 

[2 

8 1 

- TT, 

-TT 

— TT. —TT 

3 

3 J 


[3 

3 J 


(ii) For any fc € N there exist a real compactly supported scaling function ipp S 
C^(]R) and a real compactly supported associated wavelet (pM £ C^(]R) such that 
<Pf( 0) = 1 and 


/ x^(pM{x)dx = 0 

JR 


for ain G {0,1, • • • ,k — 1}. 


In both cases we have that {(fym : u £ N U 0, m € Z} is an orthonormal basis in 
L^(R) where 




ipp(t — m), ifi; = 0,TO£Z 
2^ipM{2'"~^t — m), if u £ N,m £ Z 


and the functions (pM, ‘Pf are according to (i) or (ii). 


The wavelets in the first part of the above proposition are called Meyer wavelets. 
They do not have a compact support but they are fast decaying functions and ipM 
has infinitely many moment conditions. The wavelets from the second part of the 
above proposition are called Daubechies wavelets. Here the functions ipM,PF do 
have compact support, but they only have limited smoothness. 

Definition 9.2. [16] Let s(-) £ fl 0 < q < oo and p(-) £ P(T”) with 

0 < < p"*" < OO. Denote M^- = {m : m = 0,1,2, • • • ,2^ — 1}. 

(i) Then 

= {'^Gm}ieNo,GeG3,meMj ■ ll^llf;'!(-) < oo 

p(-).g 
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where 


/ 


= 

p(-).g 


EE 

j=0 GeGi 






1/9 


Lp(-)(T")y 


(ii) For < c», we define 




.(■) < oo 


''p(-)i9(') 


where 


I p( ) = 


l/9(-) 


Vi=0 GgGJ meMj 


LP(-)(T") 


with g(-) e P’(T"). 


Theorem 9.3. Let Xi,X 2 , ■ ■ ■ be pairwise independent identically distributed ran¬ 
dom variables with E|Xi| < oo. Let = p, and Sn = Xi + X^ + ■ • ■ + X^- Then 
Sri/n pt a.s. as n ^ oo. 


The setting is to consider a random function u given by the random series 

U = (9.1) 

fc >0 

where {^fe}fc is an i.i.d. sequence and the ipk are real- or complex-valued Holder 
functions on founded open D C M" satisfying, for some a G (0,1], 

\i^k{x)-'il;k{y)\<h{a,ilJk)\x-y\°‘ x,y £ D; (9.2) 

of course if a = 1 the functions are Lipschitz. 

Theorem 9.4. Let {^fe}fc>o be countably many centered i.i.d. random variables 
with bounded moments of all orders. Moreover let {'f’k}k>o satisfy i9.S\) . Suppose 
there is some S G (0, 2) such that 

Si-='^Uk\\L’^ <oo, (9.3) 


S 2 :='^\\'tlJk\\lJh(a,ifkY < oo. (9.4) 

k>0 

Then u defined by is almost surely finite for every x G D, and u is Holder 

continuous for every Holder exponent smaller than a512. 
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